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Abstract 

The superscaling approach is applied to studies of neutral current neutrino reactions in the 
quasielastic regime. Using input from scaling analyses of electron scattering data, predictions for 
high-energy neutrino and antineutrino cross sections are given and compared with results obtained 
using the relativistic Fermi gas model. The influence of strangeness content inside the nucleons in 
the nucleus is also explored. 
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I. INTRODUCTION 



Inclusive electron scattering at intermediate to high energies from nuclei is known to 
exhibit the phenomenon of scaling and superscaling 

mum 

. At sufficiently high 

energies, typically at least 500 MeV, one sees that near the quasielastic peak the cross section 
may be analyzed in terms of a reduced response obtained by division by a suitable N- and 
Z-weighted single-nucleon electromagnetic cross section and plotted against an appropriate 
kinematic variable to see the scaling behaviour. First, when the reduced cross section is 
seen to depend only on this kinematic variable — the scaling variable — and not on the 
momentum transfer one has scaling of the first kind. Second, if the reduced cross section 
and scaling variable have been made dimensionless via removal of the momentum scale 
characteristic of a given nucleus, and the results are seen to be independent of the particular 
nuclear species, one has scaling of the second kind. When both types of scaling behaviour 
occur one says that the cross sections exhibit superscaling. In the above-cited studies the 
appropriate reduced cross sections and scaling variables have been discussed in depth. 

One finds that in the relevant energy range in the region below the quasielastic (QE) 
peak, usually called the scaling region, scaling of the second kind is found to be excellent 
and scaling of the first kind to be quite good. Above the peak scaling of the second kind 
is good; however, scaling of the first kind is clearly violated. The last occurs for well- 
understood reasons, namely, in that region one has processes other than quasi-free knockout 
of nucleons playing an important role. Specifically, the most obvious reaction mechanism is 
that of exciting a nucleon in the nucleus to a delta, which subsequently decays into a nucleon 
and a pion. Since the elementary cross section for that process is not the elastic eN cross 
section used in defining the scaling function introduced above, and since the scaling variable 
used in the usual analysis assumes the kinematics of the elastic process N — > N, rather than 
of N — > A which would now be appropriate, it is not surprising that scale-breaking occurs. 
Additionally, meson exchange current effects are known to violate the scaling behaviour, 
although from modeling in this high-energy regime jslls, IfJ, U, 3,0] their effects appear 



not to be the dominant ones. 

What was appreciated for the first time in recent work is that it is possible to pursue 
an approach where both the QE process is active (with its reduced response and scaling 
variable) and also incorporate the inelastic process in the A-region (with its corresponding 
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for various nuclear species. Of considerable 

n 

141 ] and will be discussed in the present work 



reduced response and scaling variable). We shall roughly refer to the region of excitation 
forming a peak that lies above the maximum of the QE response as the "A-peak" , although 
it should be understood that the modeling actually includes the full inelastic response on a 
nucleon (resonant plus non-resonant) for kinematics where the A(1232) is dominant. 1 In ^| 
it was shown that an excellent representation of the total inclusive electron scattering cross 
section from the scaling region up to the peak of the A-region is attained by inverting the 
procedure. Using the two scaling functions, one for QE scattering and one for the A-region, 
along with the corresponding N- and Z-weighted elastic (eN — > eN) and inelastic (eN — > 
e'A) electron scattering cross sections one finds excellent agreement with existing high- 
quality data over a wide range of kinematics anc 
importance for what was discussed in the rest of 
is the fact that the quality of the analysis requires the phenomenological scaling functions to 
be quite asymmetric, with relatively long tails extending in the direction of higher energy loss 
(positive values of the scaling variables). Such is not typically the case with most models, 
these almost always being more nearly symmetrical about the peak in the scaling function 
(see, however, [l^J where in at least one case the correct behaviour has been obtained 
in a model). This fact casts considerable doubt on most existing models for high-energy 
scattering in the QE and A regimes if high-quality results (say better than 25%) are desired. 

Having met with success in extending the scaling and superscaling analyses from the 
scaling region, through the QE peak region and into the A region, in |l4| the scaling ideas 
were inverted: given the scaling functions one can just as well multiply by the elementary 
charge- changing (CC) neutrino cross sections now to obtain the corresponding CC neutrino 
and antineutrino cross sections on nuclei for intermediate to high energies in the same region 
of excitation. Other related work is presented in jisl 3] • Given the ability of this scaling 
approach to reproduce the electron scattering cross sections, in contrast to most direct 
modeling which fails in detail to do so, we believe that such predictions for the analogous 
CC neutrino reactions should be very robust. Clearly such results are of relevance for 
on-going studies of neutrino reactions and neutrino oscillations in this intermediate-energy 
regime. 

In the present study these scaling and superscaling ideas are carried a step further to 



1 For still higher- lying excitations and DIS a different approach must be taken (see, for example, |6j). 
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include neutral-current (NC) neutrino and antineutrino scattering cross sections, in this case 
for scattering from 12 C. Specifically, the goal is to obtain results using the same analysis as 
discussed above (and in detail in 14J) for the reactions 12 C(b>,p)vX, 12 C(z/,p)z/X involving 
proton knockout and 12 C(z/, n)u\, 12 C(u, n)u\ involving neutron knockout in the QE regime, 
the A-regime being left for a subsequent study. 

A new feature emerges with such a goal in mind, however, and that arises from the 
fact that when one has an incident lepton, a scattering with exchange of a 7, W or Z°, 
and detects the scattered lepton (i.e., a charged lepton), the t-channel exchange of the 
appropriate boson is controlled. In contrast, when the scattered lepton is a neutrino or 
antineutrino, and therefore not detected, but instead a knocked-out nucleon is detected, it 

n 

is the w-channel whose kinematics are controlled (see also |17[ for discussions of this case). 
Accordingly, in the scaling analysis it is not obvious that the two types of processes are 
simply related, and therefore to apply the scaling ideas to NC neutrino and antineutrino 
scattering, in particular for discussions of differential cross sections as in the present work, 
we first have to address the issue of how the t- and u-channels are related. 

The paper is organized the following way: in Sec.|n]we begin with a basic discussion of the 
t- and w-channel kinematics involved in the semi-leptonic electroweak processes of interest 
(Subsec. Ill Aj) followed by a brief summary in Subsec. Hi Bl of the cross section formalism and 
the ideas of scaling when inter-relating t- and u-channel processes. To keep the discussions 
relatively brief in this subsection, the development of the single-nucleon NC neutrino and 
antineutrino cross sections is placed in an Appendix. For orientation in Subsec. Ill CI the 
Relativistic Fermi Gas (RFG) model is invoked and its superscaling properties summarized. 
Then in Sec. IIHI our results are presented and our conclusions are gathered in Sec. IIVI 



II. GENERAL FORMALISM FOR U-CHANNEL SCATTERING 

We begin the general discussion of how t- and w-channel semi-leptonic reactions are inter- 
related with a summary of the relevant kinematic variables in the problem. 
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FIG. 1: Kinematics for semi-leptonic nucleoli knockout reactions in the one-boson-exchange ap- 
proximation. 

A. Kinematics 

We consider general semi-leptonic quasi-free scattering from nuclei in Born approxima- 
tion. 

We start with one basic assumption that is usually presumed to be a good approximation 
in the kinematic region where quasielastic scattering is dominant, namely, that the inclusive 
cross sections are well represented by the sum of the integrated semi-inclusive proton and 
neutron emission cross sections. In doing so we are neglecting processes that occur for the 
same kinematics, but have no emitted nucleon in the final state (photon emission, deuteron 
emission, alpha emission, coherent pion production, etc., but without an emitted nucleon). 
The process of interest (see Fig. has a lepton of 4-momentum = (e, k) scattered to 
another lepton of 4-momentum K ,fl = (e', k'), exchanging a vector boson with 4-momentum 

= K^ — K'^ . The lepton energies are e = \/m 2 + k 2 and e' = \/m! 2 + k' 2 , with m (m') the 
mass of the initial (final) lepton. For NC neutrino scattering m = vnl = (assuming zero- 
mass neutrinos). Note that no assumption such as the plane- wave impulse approximation 
is being invoked at this stage. 

In the laboratory system the initial nucleus is in its ground state with 4-momentum P% = 
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(M^, 0). The final hadronic state corresponds to a nucleon (N = p or n) with 4-momentum 
P N = {En,Pn) and energy En = yj m 2 N + p 2 N plus an unobserved daughter nucleus with 
4-momentum P B = (Eb,Pb)- As usual in semi-leptonic reactions we introduce the missing 
momentum p = — p# and the excitation energy £ = Eb — E B , with E%_= \J{M B ) 2 +p 2 , 
M B being the ground-state mass of the daughter system (for details see |2|,y, |4|). 

For NC neutrino scattering we assume that the neutrino beam momentum is specified 
and the outgoing nucleon is detected. Hence pn and the angle 9 kpN (between k and Pn) 
are given. Note that the scattered lepton's 4-momentum is not specified, as would be the 
case for t-channel scattering. In analogy with the t-channel case, we can define a w-channel 
exchanged 4-momentum 

OF = K»- I* = (1) 

The above equation yields 

q = |q'| = Jk 2 + p 2 N - 2kp N cos 9 kpN ■ (2) 

For convenience in looking at the kinematics one can use a coordinate system having the 
z-axis along q', with k and pjy lying in the xz-plane. The vectors k' and p = k' — q' lie in 
a plane forming an angle with the xz-plane defined above (see Fig. |2J). 




FIG. 2: Vectors relating t-channel and ^-channel kinematic variables. 

The exclusive process illustrated in Fig. ^ is fully determined by six kinematic variables, 
which can be chosen to be (k, pn, @k PN , P, £, <$>')• The u-channel inclusive cross section 
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for (k,pN, 9k PN ) fixed is obtained by integrating over the allowed region in the (p, £)-plane 
and over the azimuthal angle, < (f)' < 2n. Again referring to Fig. El one sees that at 
fixed w-channel scattering kinematics (i.e., the triangular region bounded by k, p^r and q' 
fixed) and for a given point in the (p, £)-plane, this <$' integration corresponds to having 
the triangle bounded by p, k' and q' fixed in size and shape but rotating about the z-axis, 
namely, through the full range of the azimuthal angle . This clearly implies that the t- 
channel momentum transfer q varies and that the usual azimuthal angle <p (rotations about 
q) does not cover the full range (0, 2ir). This has consequences that are discussed in more 
detail below. 

In order to determine the integration region in the (p, £)-plane we use energy conservation, 
obtaining the following expression: 



£ = (M° A + u') 



/ m /2 + q ,2 +p 2 + 2 q'p cos 6 qlp + J {M%Y + p 2 



(3) 



Following the usual y-scaling analysis we can now examine the various curves £ = £{p) 
that result when various choices are made for cos9 g ' p . Let us call the curves £'±(p) when 



cos^'p = ±1 : 



^/m' 2 + (q'±p) 2 + J(M%f+p 2 



(4) 



From these we can proceed to find the intersections of the curves with the axis £ = 0. 
This leads to definitions for a scaling variable y' and a maximum missing momentum Y'\ 

1 



y 

y 



(M° A + uj')JA 2 x -(M° B yW 2 -q'A x 



{M A + u')JA? x -{M%yWl + q'k x 



(5) 
(6) 



where 



W 



X 



(M° A + uo') 2 - q' 2 



Ax = - 2 (w x + {M%) 2 



m 



(7) 
(8) 



Note that these are new variables and not simply related to the variables y and Y that 
ccne fr0m the familiar y - 8Calmg analy* BflflQ. The ahowed regi on „ then detained: 
for y' < one has —y' < p < Y' with < £ < £L(p), whereas for y' > one has for 
< p < y' the range S'Ap) < £ < £'-(p), and for y' < p < Y' the range < £ < £'_(p). 



When y' — one covers the largest range in missing momentum at the minimal missing 
energy and accordingly somewhere near this point the inclusive integral is expected to be 
at a maximum; namely, this kinematic point corresponds approximately to the QE peak. 

Concerning the azimuthal integration, note that kinematic variables entering the usual 
t-channel (such as the momentum and energy transfer (q, u), the lepton scattering angle B\ 
between k and k' and the solid angle defining the outgoing nucleon momentum (# gpjv , 0at)) 
all depend on cos 0' — see the above discussions of Fig. |21 Thus the integration over 
implies an integration over the azimuthal angle 0^; however, as 0' varies, the integration 
implied over 0at is not being done at constant (q, lu). Furthermore, the range over which 
the implied 07v-integration occurs is not in general the full range. This implies that the 
symmetry properties of the responses Rk cannot be used in the case of w-channel inclusive 
scattering to eliminate some of the responses (e.g., the TL and TT terms), as is the case for 
t-channel inclusive scattering. 



B. Cross sections and scaling 

Next we turn to a discussion of the basic cross sections and scaling variables involved 
in the present study. As discussed above we consider only semi-inclusive nucleon knockout 
reactions in building up the inclusive cross sections. The usual procedure is to start with 
the Plane Wave Impulse Approximation (PWIA) for the (I, I'N) cross section and integrate 
over all unconstrained kinematic variables. Final-state interactions are then presumed to 
occur after the primary electroweak interaction with a nucleon in the nucleus and so, for 
instance, a succession of (N, 2N) steps occurring during the time evolution of the high-energy 
emitted nucleon as it proceeds through the nuclear medium can cause a redistribution of 
strength in the missing-energy, missing-momentum plane (see j^fj] for recent work along these 
lines). Such processes tend to move strength from lower missing energies to higher ones and 



thereby produce an asymmetry in the scaling function, skewing it to larger values o 



energy 
also 
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loss u) or, equivalently, in the positive scaling variable direction. Other approaches 
yield an asymmetric scaling function — in agreement with experiment — when strong final- 
state interactions are incorporated, again via a shift of strength to higher missing energies. 

In contrast, in the present work where our emphasis is placed on inter- relating various 
inclusive semi-leptonic processes, and not on detailed modeling of the reaction chain, we 
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take as given the full semi- inclusive nucleon knockout cross section (i.e., given by Nature) 
and proceed to integrate to obtain inclusive cross sections. Clearly this does not imply that 
we have a full understanding of the former, only that asymptotic states may be used to 
account for all open channels and that it is not necessary to account for the entire sequence 
of steps that yields these states. 

For t-inclusive scattering, where Q^ 1 = (u, q) is constant and the final lepton is de- 
tected (as in usual inclusive electron scattering or in charge-changing neutrino reactions), 
the inclusive cross section is calculated by integrating the semi-inclusive cross section 
da / dVLk'dk 1 dVLj^dpN over the ejected nucleon (and summing over protons and neutrons), 
whereas for the w-inclusive scattering we are considering here, where Q ,tJl = (a/, q') is con- 
stant and the final nucleon is detected, one has to integrate over the final lepton. That is 
we have 

da [ dfl d ^ a (9) 

dQi-idk' J dVLk'dk'dVLjsidpN 

da = [ dn k ,dk> ^ (io) 

dVL N dj>N J dVt k idk'dVtNdpN 
for t- and u-channel reactions, respectively. These integrals can be transformed into integrals 
in the (p, £)-plane using the relations 

dVt N dp N = ( — =— ] -pdpd£d(j) N (11) 
\PnJ 1 

dtt k ,dk' = (—] -pdpdEdcp' . (12) 
\k u J q 

This leads to the following expressions for the inclusive cross sections, in the t-channel 



dVt k idk' q J Vt J J 2tt \p N J dVt k idk'dVl N dpN 
and in the w-channel 

d ° = 2 jl r pdp f d£ r d Ji (£\ <t , (14) 

dVL N dp N q' J Vu J J 2tt \k' 2 J dQ k >dk'dQ N dp N ' 

respectively. The above expressions are simply connected to one other by interchanging the 
final lepton variables with the final nucleon variables, but for the fact that the integration 
regions T> t and T> u in the (p, £)-plane are different in the two cases. The t-channel case is 

case is treated in the following section on results. 
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To this point we have made only relatively weak approximations by assuming that the 
cross sections in the quasielastic region are dominated by integrals over the semi-inclusive 
nucleon knockout cross sections. Following Q| we write the latter in terms of products of 
single-nucleon electroweak cross sections multiplied by what may be called the reduced cross 
section E: 

da 111 



;g A D v (Q*)% v w^ Q^W,Mwk,^P,£), (15) 



dVL k ,dk'dtt N dp N (2vr) 2 2e 2E~ 
where E is the energy of the struck nucleon, g is the strength of the fermion-vector boson 
coupling and Dy(Q 2 ) = (Q 2 — My) -1 is the vector boson propagator, while and w^ v are 
the usual leptonic and (single-nucleon) hadronic tensors, respectively. Clearly other sets of 
independent variables may be used as arguments of the reduced cross section (see below). 

Next we make two stronger approximations. First, we assume that the single-nucleon 
cross section varies only slowly with (p, £) and may be removed from the integrals over p 
and S. This has been verified for t-channel reactions as long as the semi-inclusive cross 
sections are peaked at low m is Sing -ener gy and miSSi „ g _t Um (see, for e,a m p,e, Q 
In particular, for the i-inclusive case the vector boson propagator can be extracted from 
the integral, and the same applies to the single-nucleon form factors appearing in w^, as 
they only depend upon Q 2 . As a consequence, in t-channel case one can verify that the 
(p, 8) dependence of the single-nucleon cross section is weak at constant (u, q) and therefore 
its mean value (namely, integrated over 0at and divided by 2n) can be removed from the 
integrations in Eq. (|13|) . The u-channel case is more complicated and is dealt with below. 

If we make this approximation we are left with 

da 



dQ^dk' 

where 



°§F{y,q), (16) 



F iV,Q) = J Pdp J ' ^T,(q,uj,6 kk ,,(j)',p,£) (17) 

depends upon the scaling variable y and the momentum transfer q QB1HQ . Note that 
the reduced cross section E occurring above would be the spectral function S(p,£), namely 
dependent only on (p,S), were the PWIA to be assumed; however, no such assumption is 
being made here. 

Second, we assume factorization in the sense that the reduced cross section appearing 
above depends only weakly on the momentum transfer q, this dependence being contained 
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mostly in the single-nucleon cross section. Note that, for instance, residual dependence in 
E on the scaling variable y is not part of the factorization assumption. Such dependence 
would not lead to any scaling violation. This means that factorization is not equivalent to 
assuming dependence only on (p, S) as in the PWIA. Clearly missing here, for instance, are 
processes involving meson-exchange currents jsj] 3] which in this sense do not factor, as 
their dependences on q are clearly not the same as those contained in the single-nucleon 
cross section which has been divided out to define the reduced cross section. However, our 
past studies of superscaling show that, for high-energy inclusive scattering at quasielastic 
kinematics, the scaling behaviour is quite well respected, with perhaps 10% or so left to 
be explained by effects such as those from MEC which should break the scaling. Indeed, 
even with relatively strong final-state interactions one finds in some modeling [l5j that the 
scaling is maintained, suggesting that the above assumption is valid, at least in the region 
of the QE peak. 

We note in passing that the violations of scaling of the first kind, namely, some residual 
dependence on the momentum transfer q, even at the level of the above equation can stem 
from two different sources: (1) the region of integration T> t depends on q and only for 
asymptotically high q does it approach a g-independent form, and (2) the reduced cross 
section £ may contain some weak dependence on q. Indeed, from the observation that the 
approach to first-kind scaling is from above, i.e. the measured reduced cross section decreases 
with q before reaching the scaling domain (see ^, for example), it appears that (2) must 
occur, and not just (1) which would imply an approach from below, since the integration 
region increases with q. 

At high energies, where the scaling idea works and scaling of first kind is reasonably good, 
we find that F(y, q) ~ F(y) = F(y, oo) and is not a function of q, in effect validating the 
factorization assumption and the quality of the approximation where a mean-value for the 
single-nucleon cross section is removed from the integrals. This was used in ^| to predict 
the charge- changing (CC) neutrino cross section: we let Nature solve for us the integral in 
Eq. fTTJ) to obtain an empirical function F(y) from electron scattering to be used in CC 
neutrino studies. 

In the li-inclusive case the above factorization is not trivial, since Q 2 varies within the 
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integration region. However, one can again assume that 



da 



F(y\ q') , 



(18) 



~ a 



where 




E ~ F(y') , 



(19) 



provided the effective NC single-nucleon cross section 




32neq' \E N 




d(j)' 
2^ 



Z MI/ (k,k>^( 



p,p N )D v (Q 2 ) 2 



(20) 



is almost independent of (p, £) for constant (k,p N , 0k PN ). This seems indeed to be the case, 
as shown from numerical studies presented below in the Results section. Then, as in jl^ . 
the empirically determined scaling function F(y') can be used to predict realistic NC cross 
sections. 

To be able to use the scaling function obtained from analyses of inclusive electron scat- 
tering data for predictions of neutrino reaction cross sections one further assumption must 
be made, namely, the domains of integration in the integrals above must be the same or at 
least very similar. In the case of CC neutrino reactions this is clearly the case except at 
very low energies for the muon case where the kinematic dependence on the muon mass is 
important in determining T> t . For NC neutrino reactions the integration domain D u differs 
to some degree from the one that enters in electron scattering, namely, T> t . In particular, 
when determining the scaling function F(y', q') with input from electron scattering which 
yields F(y, q), clearly the first step is to use the latter evaluated at y = y' and to work in the 
scaling regime where q and q' are both large enough to make the regions in the (p, £)-plane 
extend to high-p and high-£ (see the arguments for electron scattering scaling summarized, 
for instance, in j^j]). Under these circumstances the regions denoted T> t and T> u differ signifi- 
cantly only at large £ (also at large p, but there one believes the semi-inclusive cross sections 
are negligible). Accordingly, given that the semi-inclusive cross sections are dominated by 
their behaviours at low £ and low p, one expects the results of the integrations in the two 
cases, t-channel and u-channel, to be very similar, and thus the scaling functions will be 
essentially the same. Were this not to be the case, then it would be likely that first-kind 
scaling for inclusive electron scattering would not occur, in contradiction with observation. 

A further difference between the t- and ^-scattering cases should be stressed. In both 
cases the single-nucleon cross section can be expressed in terms of response functions, as 
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shown in the Appendix. However, as mentioned above, for t-inclusive processes the special 
symmetry about the q direction can be exploited to remove the TL, TT and TL' responses 
after performing the (^-integration, which simply yields a factor 2ir. In the w-channel, 
instead, the unrestricted integration over <p' yields an effective integration over which 
is not uniform and does not in general cover all of the interval (0, 2ir). As a consequence 
the TL, TT and TL' responses do contribute. As we will show later, their contribution is 
suppressed and only the TL contribution is relevant for the kinematics of interest in the 
present study. 



C. RFG and superscaling 

In this section we discuss the NC neutrino cross section in the RFG model, which corre- 
sponds to the following excitation energy 



and spectral function 



£rfg(p) = \/m 2 N + k 2 F - J m 2 N + p 2 (21) 



Srfg(p, S) = ^0(k F - P)5(E - S RFG (p)) , (22) 



where hp is the Fermi momentum and T F = \/kp + m 2 N — the Fermi kinetic energy. 

Due to the delta-function in Eq. (}2"2"|) the integration region in the (p, £)-plane simply 
reduces to a line and the lower limit in the integral over p is given by the intercept of the 
curve £rfg(p) with £'_(p) when y' < 0. When y' > it is given by the intercept of £rfg(p) 
with £' + (p) (£'_(p)) when £'±(0) < T F (£±(0) > Tp). By solving these equations it is easily 
shown that the minimum momentum required for a nucleon to participate in the reaction is 



(u) 

Vrfg 



(23) 



Pmin 

where 

A'vVy 2 + t> - k't' P '] (24) 



(u) _ rn N 
Vrfg ~ s 



r 

is the RFG y-scaling variable associated with -u-scattering (hence the index (u) to dis- 
tinguish it from the usual t-channel variable). Moreover we have introduced the dimen- 
sionless kinematic quantities k' = q'/2m^, X' = u;'/2mjv, r' = n' 2 — X' 2 and defined 
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p' = 1 — 477(1 — m f2 /m 2 N ). The sign s is 

1 



s = se;n 



r 



'T ,z p' +r - k t p 



(25) 



As in electron scattering, it is convenient to introduce a dimensionless scaling variable 

1/2 



(«) 

RFG 



m N 



7\ 



( («) \ 2 
1+ [W] _i 
\ \m N j 



(26) 



representing the minimum kinetic energy of the nucleons participating in the reaction. By 
placing the spectral function of Eq. (122(1 in Eq. (|19|l one immediately finds the RFG scaling 
function 



3k! 
TV 



d£ / ^5 (£ - £ RFG ) = \k F (l - ^l 2 G ) 9 (l - • (27) 



Providing the single-nucleon cross section is smoothly varying within the (p, £) integration 
region, the differential RFG cross section can be factorized as shown in Eq. (|18j) with the 
scaling function given by Eq. (|27j) . More realistic predictions can be given by using, instead 
of the RFG scaling function, the empirical scaling function as determined from QE electron 
scattering, as already done in Q] for charged current reactions. These are discussed in the 
next section. 



III. RESULTS 

Before presenting our predictions for the cross section, we test the validity of the scaling 
approach in the u-channel. To this end we analyze how the effective NC single-nucleon cross 
section given in Eq. (|2()jl depends on the missing momentum p and excitation energy 
£ for selected values of the kinematical variables (k, E N ,9 kpN ). To proceed, we assume 
the proton knockout case and divide evaluated in the whole (j>, £)-plane by its value 
corresponding to p — \y'\ and £ = 0. In what follows we use the cc2 off-shell prescription for 
the nucleon current and the Hohler parameterization for the single-nucleon form factors Q| , 
ignoring the strangeness content of the nucleon, unless specified otherwise. 

The results are given in Fig. El in terms of different shadings representing the regions 
where this ratio differs from unity by at most 1%, 1-2%, 2-5%, 5-10% and more than 10% 
respectively, as indicated in the top right panel. The six graphs correspond to two values 
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of the scattering angle 9k PN - 20° (top panels) and 60° (bottom panels). In each case, the 
outgoing proton kinetic energies have been selected to correspond to the regions below, above 
and close to the peak of the differential cross section. While not shown here, the results for 



neutron knockout are very similar to those for proton knockout. 
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FIG. 3: The various regions correspond to values of the ratio 1Z = (p, £) /a^ (p = \y'\,£ = 0) 
differing from 1 by at most 1% (lowest region), 1 — 2%, 2 — 5%, 5 — 10% and more than 10% 
(highest region). For this figure proton knockout has been assumed; the neutron knockout case is 
similar and not shown. For brevity, in this figure we let 6 stand for the angle 6k PN - 



The results in Fig. El illustrate the validity of the scaling approach. Only for very large 
values of the excitation energy does the effective NC single-nucleon cross section depend 
significantly on (p,S). In fact, restricting ourselves to excitation energies below twice the 
maximum value of the RFG model, S ~ 50 MeV, the dispersion presented by the ratio is at 
most ~5-10%. 

This outcome is also in accordance with the results presented in Figs. 0] (proton case) and 
El (neutron case). Here we show the neutral current neutrino (upper panels) and antineutrino 
(lower panels) double differential cross sections for scattering at 1 GeV from 12 C as a function 
of the ejected proton or neutron kinetic energy. The scattering angles have been fixed as in 
the previous figure. 
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Beginning with the RFG model, as in past work the Fermi momentum for 12 C is taken 
to be fcp = 228 MeV/c and results are given both using the full RFG model (short-dashed 
curves) and making use of the factorization approach assumed in Eq. JTHj) with the ^-channel 
NC single-nucleon cross section evaluated at p = Urfg an d £ = £rfg (solid lines). One sees 
that the two sets of results almost coincide in the whole T/v region where the RFG cross 
section is defined, indicating that the scaling argument works very well. 

Hence we may use the phenomenological scaling function extracted from (e, e') data, as 
was done in our previous CC neutrino reaction analysis Q| , to predict NC neutrino-nucleus 
scattering cross sections. These are also plotted in Figs. |U and El as long-dashed lines: they 
are seen to be lower by about 25% at the peak than the RFG results, an effect similar to 



what was found in 



14j for the charge- changing processes. Moreover, the empirical scaling 



function leads to cross sections extending both below and above the kinematical region where 
the RFG is defined. In particular, the long tail displayed for low T/y-values (corresponding 
to positive values of the scaling variable ip') is noteworthy. This tail arises not only from 
the asymmetric shape of the phenomenological scaling function, but also from the effective 
NC single-nucleon cross section, which increases significantly for low T/v values. 

On comparing Figs. 0] and 03 we see that the shapes of the cross sections for proton and 
neutron knockout are very similar, although the magnitudes are somewhat different: except 
for antineutrinos at forward angles, where the cross sections are very small, the neutron 
knockout results are 30-50% higher than for proton knockout. This occurs because (in 
absence of strangeness) both the vector and the axial-vector contributions are larger for 
neutrons than for protons, and they sum up. In particular, the AA piece is the same for p 
and n, since Ga p = —G An (see Eq. (A33)). On the other hand, from Eqs. (A28,A29,A34,A35) 
one has that Ge p — —Geu and Geu — ~Ge p - Hence, when compared with electromagnetic 
interactions, the roles of protons and neutrons are reversed in the weak neutral sector, so that 
\Geu\ » \Ge p \- Similarly, from Eqs. (A30,A31,A36,A37) one finds that |CrMn| > 
For antineutrinos things are more delicate, since the VA response has the opposite sign. For 
instance, for neutron knockout at 9 n = 20° the sum VV+AA almost exactly cancels the 
interference, explaining why the forward angle v neutron cross section is so small. 

In Fig. |U1 the contributions of the separate responses to the total RFG cross section 
are displayed. Clearly the dominant contributions arise from the Rt and Rt 1 responses, 
and in the case of neutron knockout, from R^ at low values of the kinetic energy. Note, 
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FIG. 4: (Color online) Quasielastic differential cross section for neutral current neutrino and an- 
tineutrino scattering at 1 GeV from 12 C for proton knockout obtained using the RFG (short- 
dashed), the factorized approach with the RFG scaling function (solid) and the empirical scaling 
function (long-dashed). For brevity, in this and in the following figures we let On stand for the 
angle 6 kpN ■ 



however, that while not dominant (see the discussions in the Appendix) the Rtl response 
does provide an important contribution at backward angles. In particular, since it is negative 
at low kinetic energies and positive at high, it skews the overall cross section to higher values 
of T p or T n . Such an effect is, as discussed above, absent for t-channel scattering where the 
TL, TV and TT responses are zero. 

Finally, in Figs. (proton knockout case) and |H1 (neutron knockout case) we explore the 
dependence of the cross section upon the strangeness content of the nucleon. We compare 
the results obtained from the phenomenological superscaling function in a situation where 
no strangeness is assumed (solid line) with the ones obtained including strangeness in the 
magnetic (long-dashed) and axial- vector (dotted) form factors, using for fi s = g[ s }(0) a 
representative value extracted from the recent world studies of PV electron scattering 1231 1 
and taking g s A = Gy(0) to be —0.2 



241 . The effects from inclusion of electric strangeness 
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FIG. 5: (Color online) As for the previous figure, but now showing the neutron knockout case. 



are not shown here, since G E has almost no influence on the full cross sections. 

Starting with the proton knockout results in Fig. we see that for the v case magnetic 
strangeness tends to decrease the cross section, whereas for v it has the opposite effect (the 
forward-angle v cross sections are rather small and not considered in this discussion). For 
both v and v the axial strange contribution tends to increase the cross section, and so the 
net effect of incorporating both types of strangeness content is relatively larger in the v case 
than in the v case. 

On the other hand, for the neutron knockout results shown in Fig. |H1 the situation is 
somewhat different: for v the roles of magnetic and axial strangeness are reversed from 
what is seen for proton knockout, an effect which is easily understood by examining the 
sign changes that occur in going from protons to neutrons (see Appendix). Specifically, 
Gm p and Gmu are opposite in sign, whereas, being isoscalar, is the same for protons 

(s) 

and neutrons. Similarly, being isoscalar G A does not change sign in going from protons 
to neutrons, whereas, being isovector, does. The v case is anomalous: in this case 
the interference VA response tends to cancel the VV+AA contributions. Accordingly, for 
neutron knockout including magnetic strangeness, which increases both the VV and the VA 
responses, has little net effect on the cross sections, since the two effects cancel out. 
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FIG. 6: (Color online) The separate contributions of the response functions of Eqs. ()A.15IA.17|) 
to the RFG neutrino cross section: L (short-dashed), T (solid), T' (dashed), TT (dotted), TL 
(dot-dashed), TL' (double-dashed). The upper panels are for proton knockout and the lower for 
neutron knockout. 



IV. CONCLUSIONS 



In previous work the superscaling formalism was applied to charge- changing neutrino 
and antineutrino reactions with nuclei. Using QE electron scattering data, typically at 
energies above roughly 500 MeV and up to a few GeV, those analyses resulted in a universal 
scaling function: scaling of both the first and second kinds was demonstrated for the region 
of excitation lying below the QE peak. This study was supplemented by an analysis of 
the region lying above the QE peak where the excitation of the nucleon to the A(1232) 
dominates, and again it was shown that scaling occurs in this domain, albeit with a different 
scaling function and scaling variable, as expected. Putting these together (two universal 
scaling functions and two scaling variables, together with the elementary eN — > eN and 
eN — > e'A cross sections) one has a very good representation of inclusive electron scattering 
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FIG. 7: (Color online) Effects of strangeness and radiative corrections in neutrino and antineutrino 
cross sections: no strangeness (solid), fj, s = 0.55 (dashed), g s A = —0.2 (dotted) and all the above 
effects included (dot-dashed). The case of proton knockout is assumed. 



at intermediate-to high energies from well below the QE peak up to at least the peak of 
the A-dominated region. Importantly, this high quality agreement with experiment requires 
a rather asymmetric scaling function (with a long tail extending to high energy loss) and 
from other studies undertaken by us and by others it is known that usually the results of 
modeling yield nearly symmetric scaling functions, clearly at odds with the data. It is then 
straightforward to insert, instead of the EM cross sections, the elementary CC neutrino and 
antineutrino cross sections to obtain CC cross sections on nuclei, as discussed in our previous 
work. 

In the present study we have extended that superscaling approach now to include 
quasielastic scattering via the weak neutral current of neutrinos and antineutrinos from 
nuclei at intermediate-to-high energies. The same asymmetric QE scaling function and scal- 
ing variable employed in the CC study is also used here for the NC predictions, the essential 



20 




T n [MeV] T n [MeV] 



FIG. 8: (Color online) As for the previous figure, but now for neutron knockout. 



change being simply to insert the NC neutrino and antineutrino vN and vN cross sections 
in place of the CC cross sections. Less obvious is the application of the superscaling ideas to 
the different type of inclusive reaction that must practically be considered. Namely, while 
in the CC reaction studies the relevant reaction involves an incoming lepton (a v or v) and 
detection of the corresponding charged lepton at a given scattering angle, just as in electron 
scattering with incident and scattered electrons (both are t-channel inclusive processes), the 
NC reaction is different. Here one has an incident v or i>, but now a proton or neutron 
ejected at some angle, whereas the scattered v or v is not detected — this is a so-called 
w-channel inclusive process. Thus, in the present work we have had to explore the validity of 
the superscaling ideas when applied to such n-channel scattering, again at intermediate-to- 
high energies where the scaling approach can be expected to apply. Our results indicate that 
this is the case and therefore that the scaling analysis used for CC reactions should also work 
for NC scattering. Additionally, the use of symmetries about the momentum transfer direc- 
tion in all unpolarized t-channel inclusive processes which leads to only three independent 
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response functions, L, T and T', is not applicable for w-channel inclusive scattering. There 
one has in addition the remaining responses TL, TV and TT, which cannot be eliminated 
using symmetry arguments. The results obtained in the present work show that of these 
only the TL response appears to play a significant role, at least for the kinematics chosen 
here. 

In the present work scattering of neutrinos and antineutrinos at 1 GeV from 12 C has 
been taken as representative and also since it is relevant for on-going neutrino oscillation 
experiments. Cross sections at other kinematics and for other nuclei may be obtained 
by contacting the collaboration. Several conclusions emerge from examining the results 
obtained. 

First, the NC neutrino and antineutrino cross sections are seen to be in roughly a 2:1 
ratio (v.v) at backward vN scattering angles, whereas at forward scattering angles the 
antineutrino cross sections are suppressed by an order of magnitude or more. This holds 
true for both proton and neutron knockout; moreover, the neutron knockout cross sections 
are somewhat larger than the proton knockout cross sections because of the NC single- 
nucleon form factors that enter in the two cases (see text for details). These results are 
also rather different from the corresponding CC reactions where it was observed that, for 
the kinematics chosen, the antineutrino cross sections are typically much smaller than for 
neutrinos. 

Second, the interplay of the various responses (L, T, TL, TT, T' and TL 1 ) is not trivial: 
in the various channels, v and V, proton and neutron knockout, they play different roles. For 
example, the TL response is negative at low nucleon knockout energies and positive at high 
energies, producing a shift of the total cross sections to higher energies than would occur 
with only the "usual" responses L, T and T. 

Finally, the effects of strangeness are relatively large and different for the various channels, 
implying, as in past studies, that high-quality measurements with v and v together with 
proton and neutron knockout hold the potential to yield more information on the strangeness 
content of the nucleon. 

In summary, the current study employs the same superscaling approach used previously 
for CC neutrino reactions now applying it to NC neutrino scattering in the QE region. 
Building in the correct scaling function, in contrast to direct modeling which usually fails 
to some degree when applied to electron scattering and therefore must surely fail to the 
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same degree when applied to other semi-leptonic processes, is an essential ingredient in this 
approach. In on-going work our intent is to incorporate u-channel inclusive cross sections for 
excitations in the A region and beyond as in our previous work; however, such investigations 
are more involved than the QE study presented here, since then the final state involves both 
a nucleon and a pion, and thus are postponed to the future. 
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APPENDIX: SINGLE-NUCLEON CROSS SECTION 

In this Appendix we provide the elementary cross section for the reaction 



The single-nucleon cross section a ~ l^w^ is given in term of the leptonic tensor ( assuming 
m = m' = 0) 



v{v)N 



u(v)N . 



(Al) 



(A.2) 



with x — +1 f° r neutrinos and —1 for antineutrinos, and of the hadronic tensor 



(A.3) 
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This can be decomposed into a symmetric 

< = < y + wZ (A.4) 



Kv = -ww(r) <T + + w 2V {r)X»X» (A.5) 



<a = -*>ia{t) (gr+ + w 2A (r)X^X v 



T 



-u 1A (r) + u 2A (t) { K »ri v + rfK") (A.6) 

r 

and an antisymmetric 

< = <a = 2iw 3 (T)e^ Va ^ + w 4 (r) («V - ^«") (A.7) 

tensor, where 

X » = rf + ———k^ ° n -^ hM rf + k" , (A.8) 
r 

having introduced the dimensionless variables k m = (A, re) = Q fJ, /2rriN, = P^/m^, t = 
re 2 — A 2 . Note that u XA (the pseudoscalar term), u 2A and do not contribute to l^w 1 *", 
since 

= {KpK' v + AT^ - ^„ K • K') (K** — K'^) — ( A.9) 

and 

= (KpK'j, + - <7^if ■ - = (A. 10) 

if the leptons are massless. By contracting the above tensors we get 

IfivW^ = x {v L R L + v T R T + v TT R TT + vtlRtl 

+ x(2v t >Rt> + 2v tu R tu )} , (A.ll) 

where x = 2ee' cos 2 6^/2, 9i is the lepton scattering angle, p = t/k 2 and 

10 1 

v l = P 2 , w T = -p + tan 2 — , v TT = --p (A.12) 



$Z / o 0; 1 0; / a -i j\ 

u T / = tan — y p + tan z — , v TL > = --^=ptan — . (A.14) 

The response functions are 

R L = w 00 , R T = w 11 + w 22 , R TT = w 22 - w 11 (A. 15) 

R TL = V2 (w 01 + w 10 ) (A. 16) 

R T , = i w 21 , R TL , = iy/2w 20 . (A. 17) 
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In terms of the structure functions w\, w 2 , W3 the above response functions read (for 
on-shell nucleons, r\ ■ k = r): 

K 2 

Rl = -wx{t)— +w 2 {r)(e + X) 2 (A.18) 

T 

R T = 2iwi(t) + w 2 (r)r] 2 sin 2 9 (A.19) 

R TT = -w 2 {T)r] 2 sin 2 9 cos(20) (A.20) 

-Rtl = 2V2w 2 (t)(£ + X)r) sin 6* cos (f) (A. 21) 

i? r , =2w;3(r)-(e + A) (A.22) 

Rtl> = 2^/2wz{t) kt) sin 9 cos , (A. 23) 



where the angles # and define the bound-nucleon direction with respect to the reference 
system used in the t-channel (q along the z-axis), its energy being e = E/m^. Note that in 
this system <fi = (the outgoing nucleon's azimuthal angle). 

In the usual t-channel inclusive scattering the TT, TL and TV responses vanish, since 
they are integrated over the azimuthal angle throughout the full range (0,27r); however, 
this does not occur in w-channel inclusive processes, where the integration over the outgoing 
lepton implies an integration over the full range of 0', but not of <fi. 

From Eqs. (A18-A23) we see that two of the responses are proportional to the small 
bound-nucleon momentum parameter 77 = 1/4; namely Rtl and Rtv are both 0(r)), and 
therefore vanish in the limit 77 — > 0. Accordingly, the TL and TL' responses are expected to 
be smaller than the L, T and T' responses which survive in the limit rj — > 0. On the other 
hand, since Rtt is 0(rj 2 ), one expects that the TT contributions should be the smallest, as 
is verified by examining the results in Sec. III. 

In terms of single-nucleon form factors the structure functions are (a = p,n): 

w la (r) = rG 2 Ma (r) + (1 + r)G 2 Aa (r) (A.24) 
„„ / \ G 2 Ea (r) + tG m (t) ~ 2 

w 2a {T) = — + G Aa {r) A.25 

1 + T 

w 3a (r) = G M a(r)G Aa (r) , (A.26) 
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where 
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G Ep (r) = (2-4sm 2 e w )G T E =1 (r)-4sm 2 e w G T E =0 (r)-G^(r) (A.27) 

G En (r) = -(2- Asm 2 6 w )G t e =\t) - 4 sin 2 ^G£=°(t) - G%\t) (A.28) 

G Mp (r) = (2- 4sin 2 9 w )G T A f l (r)- 4sin 2 O w G T h f {r) - G^ {r) (A.29) 

G Mn (r) = -(2 -4 sin 2 9 w )G T A f 1 (T)- 4 sm 2 9 w Glr (T)-G[ s l(r) (A.30) 

G Ap (r) = -2Gf(r) + G^(r) (A.31) 

^(t) = 2Gf(r) + G^\r) . (A.32) 

In the above 

G T E = \r) = l -[G Ep {r)+G En {r)} (A.33) 

G T E =\r) = l -[G Ep {r)-G En {r)} (A.34) 

G T m\t) = ^[G Mp (r) + G Mn (r)} (A.35) 

^ =1 (r) = \[G Mp {r)-G Mn {r)] (A.36) 

are the electromagnetic isoscalar and isovector Sachs form factors, whereas the isovector 
axial- vector form factor is given by 

Gf(r)= 1 -(D + F)G°(T) (A.37) 

with D = 1.262/1.64, F = OMD and G%(t) = (1 + 3.32r)~ 2 . 
The strangeness form factors are parameterized as follows: 

G ( E \r) = Ps tG°(t) (A.38) 

GS?(r) = VsG°(t) (A.39) 

4 S V) = fliG5(r) , (A.40) 

with G${t) = (l + 4.97r)- 2 . 
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